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Using Monte Carlo simulation methods in the canonical and grand canonical ensembles, we discuss the melting 
and the formation of ordered structures of mixed Ar-Xe submonolayer films on graphite. The calculations have 
been performed using two- as well as three-dimensional models of the systems studied. It is demonstrated 
that out-of plane motion does not affect the properties of the adsorbed films as long as the total density is 
not close to the monolayer completion. On the other hand, close to the monolayer completion, the promotion 
of particles to the second layer considerably affects the properties of mixed films. It has been shown that the 
mixture exhibits complete mixing in the liquid phase and freezes into solid phases of the structure depending 
upon the film composition. For submonolayer densities, the melting temperature exhibits non-monotonous 
changes with the film composition. In particular, the melting temperature initially increases when the xenon 
concentration increases up to about 20%, then it decreases and reaches minimum for the xenon concentration 
of about 40%. For still higher xenon concentrations, the melting point gradually increases to the temperature 
correspondingto pure xenon film. It has been also demonstrated that the topology of phase diagrams of mixed 
films is sensitive to the composition of adsorbed layers. 
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1. Introduction 

It is now well known that monolayer films of both argon and xenon on graphite form incom- 
mensurate two-dimensional solid phases at low temperatures. However, the melting transition of these 
two solid phases has a different mechanism |9]. Experimental data Ild-[l3l as well as computer simula- 
tions jlilES have demonstrated that submonolayer films of argon melt via continuous phase transition 
at the triple point temperature equal to Tt =; 49.7 K |17J- At the temperatures well below the melting point, 
the solid-like argon submonolayer fUm is rotated by about 2-3 degrees with respect to the i?30° axis of 
the commensurate ( \/3 x \/3)R30° structure I161 Ii8h20|i . Since the argon atoms are rather small, the solid 
phase is compressed with respect to the (\/3 x \/3)R30° commensurate structure, and the monolayer 
density is about 25% higher than the density of a perfect commensurate structure. 

The melting of incommensurate submonolayer fUms of xenon is of first order (Uli^l, and occurs at 
the triple point temperature of about 100 K. However, the xenon atoms are larger so that the incommen- 
surate solid phase is dilated with respect to the commensurate ( \/3 x \/3)i?30° structure. 

The differences in the structure of low temperature phases and the phase behavior of argon and 
xenon submonolayer films on graphite lead to a quite complex phase behavior of mixed films. The X- 
ray scattering studies of the Ar-Xe mixture adsorbed on graphite |23- 25] have demonstrated that the 
structure of submonolayer and monolayer solid-like phases strongly depends upon the film composition. 
Three different solid-like phases have been found. Apart from the compressed argon-like and dilated Xe- 
like incommensurate phases, being stable for sufficiently low and high xenon mole fraction, respectively, 
the formation of krypton-like ( \/3 x \/3)R30° commensurate structure has been found over a rather 
wide range of the mbcture composition. The calorimetric study of Ma et al. | 26] has shown that even very 
small amounts of xenon, about 1.5%, added into submonolayer argon films, lead to the disappearance of 
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the heat capacity peak attributed to the orientational transition in the argon-like incommensurate solid 
phase |16, 27]. 

A vast majority of theoretical studies of mixed adsorbed layers has been based on lattice gas models 
(i^ - fiill . which do not constitute a good basis for the discussion of the incommensurate-commensurate 
transitions in the films of rare gases on graphite. Such models cannot properly describe the incommen- 
surate floating solid. However, there have also been some attempts to construct theoretical models for 
the commensurate-incommensurate transitions in monolayer films of rare gas mixtures on graphite 
[36.1 . The primary aim of the models proposed by Marti et al. |34, 36], was to explain the anomaly in the 
phase behavior of Ar-Xe and Kr-Xe submonolayer films on graphite. Namely, experiments have demon- 
strated fi^l that less krypton than argon is needed to induce the formation of commensurate phase. On 
the other hand, a rather general mean-field model of Villain and Moreira [35] requires the introduction of 
several approximations in order to make the resulting equations numerically tractable. Consequently, the 
agreement with experimental data is rather poor. Nevertheless, these authors have derived qualitative 
phase diagrams for the Ar-Xe mixture adsorbed on graphite. The theory predicts the existence of incom- 
mensurate, Ar-hke and Xe-like structures at low temperatures, and the formation of the commensurate 
structure at higher temperatures over a rather limited range of the Xe mole fraction. 

In this work, we present and discuss the results of rather extensive Monte Carlo simulations of mixed 
Ar-Xe films on graphite. Our main goal has been to investigate the structure of low temperature solid 
phases as well as to determine the changes of the melting temperature with the mixture composition. 
However, we also discuss the evolution of phase diagrams resulting from the changes in the film compo- 
sition. 

The paper is organized as follows. In the next section we present the model used and describe the 
Monte Carlo method used to determine the properties of mixed submonolayer films. Then, in section 3 
we briefly discuss the behavior of pure Ar and Xe films. The last section 4 is devoted to the presentation 
of the results for the mixed Ar-Xe films. 



2. The model and Monte Carlo methods 

The interaction between adsorbate atoms is assumed to be represented by the (12,6) Lennard-Jones 
potential 

Hjinj)-48ij \[atjlrijY^-{ai,j/njf], (2.1) 

where rij is the distance between a pair of atoms and i and j mark the species Ar and Xe. The values 
of the parameters and atj used in this work are given in table [l] The corresponding parameters 

Table 1. Lennard-Jones parameters for Ar and Xe used in this work. 
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K 


Ar,Ar 


3.4 


120.0 


Kr,Kr 


4.1 


221.0 


Ar,Kr 


3.S 


162.85 



representing the Ar-Xe interaction, also given in table [T] have been obtained using the usual Lorentz- 
Bertholot combining rules: 

1 , , , 

cTi.j = - [cri.i + o-;,;J and Eij = y/eij^ej^. (2.2) 

The potential i2.1\ has been cut at the distance Satj. 

We are aware of some drawbacks that the assumption of the LJ potential has got, and that other 
authors have used different potentials to reproduce experimental data for adsorbed films of pure rare 
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gases on graphite fsT*]. Also, we have not taken into account the surface mediated interactions, which are 
known to affect the strength of adsorbate-adsorbate interaction in the vicinity of solid substrates |7J. 

The interaction of rare gas atoms with the graphite basal plane can be represented by the potential 
proposed by Steele (H] 



Vi{x,y,z] 



-gs,i 



[ Vo,i (z) + ^ Vkj{z)ft{x,y]j, i = Ar or Kr. 



(2.3) 



In the above, the first term in the square brackets is the fluid-solid potential averaged over the entire 
surface, while the second term represents the corrugation part of the fluid-solid potential. Assuming 
that the interaction between an adsorbate atom and the carbon atom of the graphite substrate is also 
represented by the Lennard-Jones potential, the explicit expressions for VQ iiz), the Fourier components 
Vjc^iz) and the functions fkix,y) are given by the following equations: 



47i:A° °5 
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and 



|(|)V3M-2(|)r,M 



fkix,y) = ^exp[iqfc,,T], t = (x,y) 



(2.4) 



(2.5) 



(2.6) 



with the sum running over all graphite reciprocal lattice vectors of the length q^. In the above equations 
- o'i,c/'2i> where a\ - 2.46 A is the graphite lattice constant, i.e., the distance between the centers 
of adjacent carbon hexagons, the values of at e and Sggj (i - Ar, Xe) are given in table [2] Az - 3.4 A is 
the spacing between graphite planes, - 5.24 A^ is the area of the graphite unit cell, K2 and Ks are the 
modified Bessel functions of the second kind and of the second and fifth order respectively, and q^'s are 
the lengths of the graphite basal plane reciprocal lattice vectors. 



Table 2. The parameters describing the Ar-graphite and Xe-graphite interaction, obtained using the 
Lorentz-Bertholot combining rules [given by equation (2.2H and assuming that ^ = 28 K and aQ Q = 
3.4 A. 



/ 




Ai 




K 


A 


Ar 


58.00 


3.40 


Kr 


78.66 


3.75 



In the case of only partially filled monolayer films and at sufficiently low temperatures, the promo- 
tion of the second layer is likely to be negligibly small. This allows us to consider a simple strictly two- 
dimensional model with the external field of the form 

v[x,y) = V\,jfi{x,y) = -Vb,; {cos(qir) -1- cos(q2r) -1- cos([qi - q2]r)} , (2.7) 

where the parameter 14),; (' - Ar, Kr) determines the amplitude of the corrugation potential. The mag- 
nitudes of 14), i (Vb,Ar - 0-07 and Vti,xe - 0.08) have been adjusted in such a way that the results for each 
component are more or less consistent with the full 3D calculations. 

Simulations have been performed using the Monte Carlo method in the canonical and grand canonical 
ensembles [39, 40]. In the case of two-dimensional model, the rectangular simulation cell of the size 
Lai ^ L v^fli 12, with L = 60 and with the standard periodic boundary conditions has been used. In three 
dimensional calculations, the simulation cell has been the rectangular parallelepiped of the size 60fli x 
60fli v^/2 X lOfli, with the periodic boundary conditions applied in the directions parallel to the substrate 
surface and with the reflecting hard wall located at z = IOai . 
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The quantities recorded included the average potential energy, {e), the contributions to the potential 
energy due to the fluid-fluid interaction, (egg) and the contributions due to the fluid-solid interaction for 
each component <egs,i> and the heat capacity obtained from the fluctuations of the potential energy. 



(2.8) 



In order to monitor the structure of solid phases we have used radial distribution functions, gij{r), 
for different pairs of species i and j, and appropriate order parameters. The formation of hexagonally 
ordered phases has been monitored using the bond-orientational order parameters (41II42I1 



1 



(2.9) 



measured separately for each adsorbate (i - Ar or Xe). In the above, the first sum runs over all atoms of 
the /-th component, the second sum runs over all nearest neighbors of the same type, (pm,n is the angle 
between the bond joining the atoms m and n and an arbitrary reference axis, chosen here to be the x- 
axis of the simulation cell, and iVb,; is the number of bonds between pairs of the Uke atoms. Also, we have 
calculated the total bond-orientational order parameter 



TrELexp(i6(^'m,n) 

^^b m n 



(2.10) 



where the first sum runs over all atoms in the system and the second over all nearest neighbors. 

The above defined bond-orientational order parameters make it possible to detect the hexagonally 
ordered structures, but are not suitable to distinguish between the commensurate and incommensurate 
phases. In the commensurate phase, the atoms are localized over the centers of carbon hexagons, and the 
appropriate order parameter allowing to monitor such localized structures can be defined as Ii43il 



1 ^ 

oJ^i m n=l 



(2.11) 



The first sum is taken over all atoms of the i-th component, while the second sum runs over the six recip- 
rocal lattice vectors q„ of the graphite substrate and rm,/ is the position of the m-th atom of component /. 
The above defined order parameters have been supplemented by the corresponding susceptibilities 



Iop=^[<Op2>-<op>2], 



kT 



(2.12) 



where 'op' stands for any of the above given order parameters. 

When grand canonical simulations have been carried out, we have also recorded the adsorption- 
desorption isotherms. 

Throughout this paper, we use reduced quantities, assuming that the graphite lattice constant ai 
is the unit of length, and the Lennard-Jones parameter eAr,Ar is the unit of energy. We have decided, 
however, to give the temperature in Kelvins, as it allows for easier comparison of our simulation results 
with experimental data. All the densities are expressed in commensurate monolayers. 



3. The results for pure Ar and Xe films 

The phase diagram of xenon monolayer on graphite exhibits the vapor-liquid critical point and the 
triple point (see parts (a) and (b) of figure [!}■ The triple point temperature, equal to Ttr ~ 89 K, agrees 
quite well with some experimental data |44], yet it is lower than the value of about 99 K stemming from 
other experiments |22, 45] and from the recent Monte Carlo results of Przydrozny and Kuchta |37]. The 
underestimation of the triple point temperature is associated with our choice of the interaction poten- 
tial and its parameters. Przydrozny and Kuchta applied a semi-empirical potential proposed by Aziz 
and Slaman i46ll . while we have used a simple Lennard-Jones potential. We should mention that in the 
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Figure 1. The phase diagrams for pure Xe (parts (a) and (b)) and pure Ar (parts (c) and (d)) monolayer 
films on graphite derived from grand canonical Monte Carlo simulation. Parts (a) and (c) show the 
temperature-chemical potential projections, while parts (b) and (d) show the temperature-density pro- 
jections. 



earlier molecular djmamics simulation studies of the melting transition of xenon on graphite (i^, fi^l . 
also based on the Lennard-Jones potential, but with slightly different values of the parameters £xe,xe 
and CTxe. the triple point temperature was found to be located at Ttr ~ 0.4£xe,xe/fc. while our result is 
Ttr ~ 0.402£xe.xe Ik- Also, the critical point temperature Ta ~ 109 K is lower than the experimental value 
by about 18 K |45i. Nevertheless, the qualitative agreement with the available experimental data is good 
enough to assume that the results for the mixed films are also qualitatively correct. 

In the case of argon, the phase diagram derived from our grand canonical simulation (see parts (c) 
and (d) of figure[l) agrees very well with experiment. In particular, the triple point temperature, equal to 
49.5 + 0.5 K, is practically the same as the experimental value of 49.7 K [111 . Also, the critical temperature 
agrees very well with experiment ligll . 

The freezing of submonolayer xenon and argon films leads to the formation of incommensurate struc- 
tures. The xenon incommensurate solid attains the density of about 0.85 at the monolayer completion and 
is expended with respect to the commensurate ( ^3 x \/3}R30 structure. Figure|2]shows the temperature 
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Figure 2. The temperature changes of the order parameters ^6,Xe ^nd Oxe for Pure Xe films on graphite 
at different densities (given in the figure). The results for pc = 0.5 and 0.7 have been obtained using a 
two-dimensional model, while those for pc = 0.8 and 0.9 using a three-dimensional model. Filled and 
open circles mark the bond-orientational and positional order parameters, respectively. 



23601-5 



A. Patrykiejew 



changes of the order parameters ^6,xe and (Pxe for xenon films of different total densities and one sees 
that the order parameter Oxe is quite low even at very low temperatures, indicating the lack of localiza- 
tion of adatoms over the minima of the graphite lattice. On the other hand, the bond-orientational order 
parameter demonstrates the formation of hexagonally ordered phase below the freezing point. 




Figure 3. The snapshot of configuration for pure argon film recorded for pc = 1.0 at T = 10 K. The dashed 
lines show that the film is rotated with respect to the symmetry axis of the commensurate phase by the 
angle a = 3°. 

The argon incommensurate solid also shows a well developed hexagonal symmetry, but it is con- 
tracted with respect to the commensurate ( ^3 x v^)i?30 structure and attains the density of about 1.25. 
At sufficiently low temperatures, the film exhibits epitaxial rotation of about 3 degrees (see figure [3). 
Again, this result agrees very well with earlier theoretical and computer simulation IigIIztI] re- 

sults. 

The results of Monte Carlo simulation for pure argon and xenon films given above will serve as ref- 
erence data for the study of mixed films. 



4. The results for mixed films 

We begin with the presentation of canonical ensemble Monte Carlo simulation results aiming at the 
determination of the melting temperature and the structure of solid phases in submonolayer mixed films. 
Since the solid phases (commensurate and incommensurate) exhibit hexagonal symmetry, the location of 
the melting point can be estimated using the bond-orientational order parameter, I'e. and its susceptibil- 
ity, x^ti ■ Of course, one also expects that the melting transition is manifested by sudden changes of the po- 
tential energy and the heat capacity anomalies. Figure |4]gives an example of our results, obtained for the 
submonolayer film of the total density pc - 0.4 and the xenon mole fraction equal to xxe - 0.1667. Part (a) 
of figure|4]shows the heat capacity curve and one sees a sharp peak at the melting point at T =; 54 K. At 
the same temperature, the total potential energy u and the contributions to the potential energy due to 
Ar-graphite and Xe-graphite interactions exhibit sudden drops (see part (b) of figure |4). Finally, part (c) 
of figure |4] which shows the temperature changes of the bond-orientational order parameter and 
its susceptibility demonstrates that the melting transition is accompanied by the loss of hexagonal 
ordering. It should be emphasized that a large increase of the Ar-graphite and Xe-graphite interaction 
energies accompanying the freezing transition marks a sudden increase of localization of the adsorbed 
argon and xenon in the solid phase. Upon a decrease of temperature, the localization of xenon gradually 
increases, whUe argon exhibits a decrease of localization at temperatures below T =; 30 K. This behavior 
can be attributed to the transition between the commensurate phase, stable at T > 30 K, and the incom- 
mensurate phase, stable at r < 30 K. Note that the transition is not accompanied by any changes in the 
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behavior of the bond-orientational order parameter, 
but produces a well seen heat capacity anomaly. 

The inspection snapshots of configurations recorded 
during the simulation runs have shown that the 
commensurate phase is mixed, while the incommensu- 
rate phase consists of argon only. In the case of smaU 
xenon mole fraction, as in the system considered now, 
we expect to observe only a partially developed com- 
mensurate phase. Indeed, the snapshot given in fig- 
ure [5] (a) shows that the film peripheries are predom- 
inantly occupied by argon atoms, which also show a 
rather high degree of incommensuration. At the tem- 
perature below commensurate-incommensurate tran- 
sition, we find coexisting domains of mixed commen- 
surate and argon-like incommensurate phases (see fig- 
ure[5](b)). In the snapshots given in figure [s] we have 
assigned the atoms to commensurate and incommensu- 
rate positions using the following order parameter (s^ : 

(/)(r) - cos(qir) -1- cos(q2r) -1- cos([qi - q2]r) , (4.1) 



40.0 

T[K] 



and assuming that the atom is commensurate (incom- 
mensurate) when > ((/) ^ O).0ne should note that 
even in a rather small system used, consisting of only 
480 atoms, the argon-like incommensurate domain ex- 
hibits epitaxial rotation, just the same as observed for 
pure argon films. 

In order to determine the locations of the commen- 
surate-incommensurate transition, we have monitored 
the behavior of the order parameters Oat and <I>xe, defined by the equation <2.11) . Figure [6] shows the 
changes of these two order parameters with the xenon mole fraction at two different temperatures in 



Figure 4. The temperature changes of the heat 
capacity (part (a)), the total potential energy 
and contributions to the potential energy due 
to for the mixed submonolayer film of pc = 
0.4 and xxe = 0.1667, obtained using two- 
dimensional model. 





Figure 5. The snapshots obtained for the film of the total density pc = 0.4 and xxe = 0.1667 at T = 36 (left 
panel) and 24 K (right panel). Black dots mark the centers of graphite cells, open circles with thin and 
thick lines represent argon atoms being commensurate and incommensurate with the graphite lattice, 
whUe larger light shaded and dark shaded circles are the xenon atoms being commensurate and incom- 
mensurate with the graphite lattice. The dashed line in part (b) shows that the argon-like incommensurate 
phase exhibits epitaxial rotation. 
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Figure 6. The order parameters c&Ar and Oxe versus the xenon mole fraction in submonolayer film of the 
density pc = 0.4 at two different temperatures, shown in the figure. The dotted vertical lines mark the 
regions of xxe over which the pure mixed commensurate phase appears, and the vertical dash-dotted 
line marks the xenon mole fraction at which the domain of commensurate phase disappears. 



submonolayer films of the total density pc - 0.4. Quite similar results have been obtained for the films of 
different total densities and using two- as well as three-dimensional models. From the observed changes 
of the order parameters <l)Ar and <l>xe> it follows that the increase of the xenon concentration leads to a se- 
quence of changes in the film structure. For small xxe we find that xenon is highly localized, while the de- 
gree of localization of argon increases with jcxe- In this region, the film consists of two coexisting phases: 
one being the incommensurate argon-like solid and the second being the mixed commensurate solid. A 
gradual increase of the xenon mole fraction causes the commensurate domain to become larger and the 
size of incommensurate domain to shrink gradually. Then, there is a region of xenon concentration over 
which both adsorbates are highly localized. This corresponds to the presence of pure mixed commensu- 
rate phase and terminates at the xenon mole fraction close to about 0.4. Then, both order parameters 
gradually decrease when xxe increases up to about 0.75. In this region, the mixed commensurate phase 
coexists with the demixed xenon-like incommensurate phase. Finally, for xxe exceeding about 0.75, the 
film consists of xenon-like incommensurate phase with the argon atoms located at its peripheries. This 
has been confirmed by the inspection of snapshots and radial distribution functions. 

The central result of the canonical ensemble Monte Carlo study is given in figure[7] which contains the 
phase diagrams showing the locations of the melting transition and the regions of stability of different 
solid phases in the fUms of different total densities. Part (a) of figure [7]gives the results for submonolayer 
films of different total densities, equal to 0.4, 0.667 and 0.8. The results for pc - 0.4 and 0.667 have been 
obtained using a two-dimensional model, while those for pc = 0.8 have been obtained within a more re- 
alistic three-dimensional model. The locations of the melting point are more or less the same over a wide 
range of xxe between and about 0.8. The independence of the melting temperature of the total density 
indicates that the melting occurs at the triple point temperature. For the xenon concentration higher than 
0.8 the triple point melting occurs for pc - 0.4 and 0.667, but not for pc - 0.8. This suggests that the den- 
sity Pc - 0.8 is higher than the liquid density at the triple point. A rather sharp increase of the melting 
temperature with the xenon mole fraction, for xxe above 0.8, results from the fact that the increase of 
xenon concentration brings the film closer to the monolayer completion. We should emphasize that even 
for xxe close to unity there is no trace of the promotion of adsorbed argon and xenon to the second layer, 
even at the temperatures above the melting point. Thus, the films remain practically two-dimensional. 
Although we have not performed any simulation at pc - 0.8 using the two-dimensional model, it can 
be anticipated that the results should be quite the same as those obtained with the three-dimensional 
model. The two-dimensional approximation is expected to faU when the film density starts to exceed the 
monolayer capacity. 
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Figure 7. The phase diagrams derived from the canonical ensemble simulations for submonolayer 
(part (a)) and monolayer (part (b)) mixed films of argon and xenon. In part (a), the filled, shaded and 
open circles show the melting temperatures for the films of different total density equal to 0.4, 0.5 and 
0.8, respectively. The results for pc = 0.4 and 0.5 have been obtained using two-dimensional model, while 
those for pc = 0.8 have been obtained using three-dimensional model. The filled squares, diamonds and 
triangles mark the stability regions of differently ordered solid phases. In part (b) the filled and open 
circles mark the melting points obtained using three- and two-dimensional models, respectively. Open 
squares give the onset of the second layer promotion and the filled diamonds show the limit of stability 
of the mixed incommensurate phase. Dotted lines show the approximate locations of stability limits of 
different solid phases. 



We have carried out the canonical ensemble calculations assuming that pc = 10. This density is lower 
than monolayer capacity of pure argon film, but it is weU above the monolayer capacity of pure xenon 
film. Part (b) of figure[7]shows the phase diagram obtained. In this case, the two-dimensional model works 
weU only in the region of the xenon mole fraction not greater than about 0.2, and starts to overestimate 
the stability of the solid phase for higher xenon concentrations. Three-dimensional calculations have 
shown that argon is partially promoted to the second layer when the temperature becomes high enough. 
The temperature at which the second layer promotion begins depends upon the fUm composition and it 
is higher than the melting temperature only for xxe not exceeding about 0.2 and lies below the melting 
temperature for higher xenon mole fractions. It is therefore not surprising that the two-dimensional 
model works well only for xxe ^ 0.2. One of the consequences of the promotion of argon atoms to the 
second layer is the increase of xenon concentration in the first layer with respect to the nominal xenon 
concentration in the simulation cell. The locations of phase transitions in figure [7] (b) have been plotted 
for the values of Xxe corresponding to the actual xenon concentration in the first layer. 

One sees that the xenon mole fraction range over which the commensurate phase is stable is con- 
siderably wider than in the previously discussed fUms of lower total density (cf. figure[7](a)). Moreover, 
we find that over a certain range of xenon concentrations, between about 0.48 and 0.83, the mbced in- 
commensurate phase appears at the temperatures just below the freezing transition, whereas no trace 
of such a phase has been found in submonolayer fUms. Upon the decrease of temperature, this phase 
transforms either into the a commensurate phase, when xxe is lower than 0.6, or into the coexisting 
commensurate and xenon-like incommensurate phase, when xxe is higher than 0.6. Figure [8] shows the 
changes of the heat capacity (part (a)) and the order parameters (part (b)) in the case of the film with 
xxe - 0.7. The heat capacity has two pronounced anomalies. The first one, at T =s 108.5 K, is the signature 
of freezing transition, accompanied by the development of hexagonal order in the film (see the behavior 
of the bond-orientational order parameter in figure|8](b)). At the temperatures between the freezing 
point and the second heat capacity anomaly at T =; 70 K, the solid phase is incommensurate and the order 
parameters <1>ai and (Pxe remain very small, as expected for the incommensurate phase. The inspection 
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of snapshots has shown that the incommensurate 
phase is mixed and the argon is partially pro- 
moted to the second layer. The second heat ca- 
pacity anomaly, at T =; 70 K, is due to the onset 
of the transition accompanied by a rather large 
increase of the order parameters Oai and <l>xe) 
indicating the increase of localization of Ar and 
Xe upon the lowering of temperature. Figure [9] 
shows the snapshots recorded at 72 K and 12 K, 
which demonstrate that the transition observed 
leads to the formation of domains consisting of 
the mixed commensurate and demixed Xe-like in- 
commensurate phases. The snapshot recorded at 
12 K also demonstrates that argon is partially pro- 
moted to the second layer, and forms a compact 
island of a solid-like phase. It is also noteworthy 
that the solid-like patch of argon in the second 
layer is located over the demixed xenon domain 
rather than over the domain formed by the mixed 
commensurate phase. This can be readUy under- 
stood by taking into account the magnitudes of 
Ar-Ar and Ar-Xe interaction energies, measured 
by the Lennard-Jones potential parameters CAr.Ar 
and eAr,xe. and of course eAr.xe is considerably 
larger than £Ar,Ar (cf. table [1). When the argon 
atoms from the second layer are located over the 
pure xenon patch, each of them has three xenon atoms from the first layer as nearest neighbors. On 
the other hand, if the argon patch were located over the mixed commensurate patch then some nearest 
neighbors from the first layer would be argon atoms, and this situation is energetically less favorable. 

When the xenon concentration exceeds about 0.83, the first layer consists only of xenon, while all 
argon atoms are promoted to the second layer. Of course, when the amount of xenon in the film exceeds 
the monolayer capacity of pure xenon fUm, then the excess of xenon is also located in the second layer. 




Figure 8. The temperature changes of the heat ca- 
pacity (part (a)) and of different order parame- 
ters (part (b)) obtained for the system with pc = 
1.0 and xxe = 0.7 using three-dimensional model 
The vertical dashed lines mark the freezing and 
commensurate-incommensurate transitions. 




Figure 9. The snapshots for the mixed film of the total density pc = 1.0 and xxe = 0.7 at T = 72 K (left 
panel) and 12 K (right panel). Black dots mark the centers of graphite cells, open circles represent argon 
atoms being commensurate with the graphite lattice, while larger light shaded and dark shaded circles 
are the xenon atoms being commensurate and incommensurate with the graphite lattice. Filled circles 
stand for argon atoms located in the second layer. 
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Figure 10. The examples of adsorption isotfierms at r = 60 K, obtained for /i^^ = -19.0 (part (a)) and 
-18.6 (part (b)). Circles show the total adsorption, while squares and diamonds denote the adsorption of 
argon and xenon, respectively. 



We now proceed to the discussion of the changes in the phase behavior resulting from grand canon- 
ical Monte Carlo simulation. The calculations have been performed under the condition of the fixed 
chemical potential of xenon, so that the xenon concentration in the film was not conserved. Along the 

adsorption isotherms obtained by changing the 
chemical potential of argon, the amounts of xenon 
change as well. Figure [lO] shows two examples of 
adsorption isotherms, both recorded at T = 60 K, 
but with different values of the xenon chemical 
potential. It is quite evident that the gas-liquid 
transition is accompanied by a sudden increase 
of the xenon density and that this effect is much 
stronger when the chemical potential of xenon is 
higher. We have determined the phase diagrams 
for a series of systems with different values of 
the xenon chemical potential, /ixe - -20.0, -19.5, 
-19.0 and -18.6. 

In the case of /ixe - -20.0, the amounts of 
xenon in the film are very low, with xxe < 0.01, 
over the entire range of temperatures and fUm 
densities studied. It is, therefore, not surprising 
that the phase diagram obtained is very similar 
to that of pure argon film (see figure ITT). In par- 
ticular, the melting transition appears to be con- 
tinuous and the solid phase is an incommensurate 
argon-like phase. However, we find that even very 
small amounts of xenon shift the locations of the 
triple and critical points towards higher tempera- 
tures. 
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Figure 11. The phase diagram for the system with 
/iXe - -20.0 (fiUed symbols) and of pure argon film 
(open symbols). Parts (a) and (b) show the temper- 
ature - argon chemical potential and the tempera- 
ture - total density projections, respectively. 
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Figure 1 2. The phase diagrams for the systems with ;Uxe = - 19.5 (parts (a) and (b)) and - 19.0 (parts (c) and 
(d)). Parts (a) and (c) show the temperature - argon chemical potential projections and parts (b) and (d) 
the temperature - total density projections, respectively. In parts (b) and (d), fiUed circles show the total 
density and open circles show the argon density. The vertical dash-dotted lines mark the temperature 
above which the commensurate solid looses stability. 



An increase of [ixe to -19.5 leads to some qualitative changes in the phase diagram topology (see 
parts (a) and (b) of figure [12). In particular, at the temperatures below about 50.5 K, the two-dimensional 
gas condenses into the commensurate Icrypton-like structure, which undergoes a transition into the in- 
commensurate phase when the argon chemical potential becomes high enough. This commensurate- 
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Figure 13. The argon-argon radial distribution functions for the system with ^xe = -19.0, recorded at 
r = 51 K (part (a)) and T = 57 (part (b)) for the values of the argon chemical potential given in the figure. 
The vertical dotted lines marlc the locations of subsequent maxima in a perfectly ordered commensurate 
phase. Part (a) also shows the argon-xenon radial distribution function obtained for /i^r = -9.9. 
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incommensurate transition is continuous. One should note that the xenon concentration in the commen- 
surate phase is rather low (xxe < 0.2), and becomes still lower in the incommensurate phase. For still 
higher value of /ixe equal to -19.0 the phase diagram topology remains the same and only the stability 
region of the commensurate phase becomes wider and extents up to T =j 57 K. Also, the commensurate- 
incommensurate transition occurs at higher values of the argon chemical potential. At the temperature 
between the upper limit of the commensurate phase, i.e., r =; 57 K, and the critical point, the gas phase 
condenses into the liquid phase. The picture presented above is very well confirmed by the behavior of 
radial distribution functions. Part (a) of figure [131 gives the argon-argon distribution functions recorded 
at the temperature of 51 K and for the argon chemical potentials below and above the commensurate- 
incommensurate transition. It is evident that the maxima, apart from the first one, coincide very nicely 
with the locations of subsequent neighbors in the commensurate phase. The stability of commensurate 
phase is due to the presence of argon-xenon nearest neighbors, and the argon-xenon distribution func- 
tion (also shown in figure[T3las a dashed line) exhibits the first maximum very close to \/3ai, as expected 
for the commensurate phase. 

At the temperature of 38 K, i.e, above the upper limit of the commensurate phase stability, the argon- 
argon distribution function recorded at /xai - -9 demonstrates the presence of a liquid phase, while at 
/iAr - -8.8 it corresponds to an incommensurate solid phase. The liquid is of course partially ordered due 

One should also note a gradual increase of 
the critical temperature resulting from the in- 
crease of xenon concentration. The phase behav- 
ior changes when the chemical potential of xenon 
is increased to -18.6. Figure [141 shows that the 
gas condenses into a liquid phase of rather high 
xenon concentration, ranging between xxe ~ 0.67 
at r = 56 K and xxe ~ 0.44 at T = 66 K. When the 
argon chemical potential increases, we observe 
the transition between the liquid and commensu- 
rate phases. This transition, quite well illustrated 
by the change in the behavior of the argon-argon 
radial distribution function given in figure [15] oc- 
curs only at the temperatures lower than 61 K. 
A further increase of the argon chemical poten- 
tial does not lead to the transition between the 
commensurate and incommensurate solid phases, 
as in the previously considered cases, but rather 
again to the liquid phase. The liquid undergoes 
a transition into the incommensurate solid-like 
phase at stiU higher values of the argon chemi- 
cal potential (cf. figure [15). This re-entrant behav- 
ior can be understood by taking into account that 
the upper limit of the film density in the commen- 
surate phase is equal to 1.0, while the transition 
into the incommensurate solid in the argon rich 
film and at the temperatures used occurs at the 
densities well above unity. An increasing chemi- 
cal potential leads to a gradual removal of xenon, 
so that the dense film becomes more and more 
argon-like. Note that the liquid-incommensurate solid transition in pure argon film at T = 56 K occurs 
at the density of about 1.07 (cf. figure [l](b)), and at still higher densities at higher temperatures. 

Concluding, we would like to present the comparison of our results with the available experimen- 
tal data for submonolayer films of the total density equal to pc - 0.4. One readily notes a qualitative 
agreement between Monte Carlo and experimental results. However, the present simulation predicts a 
considerably narrower range of xenon concentrations corresponding to the stability region of the com- 



to the effects of periodic corrugation potential. 
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Figure 14. The phase diagram for the system with 
/iXe = -18.6. Parts (a) and (b) show the tempera- 
ture - argon chemical potential and the tempera- 
ture - total density projections, respectively. The ab- 
breviations G, L, C and IC stand for the gas, liquid, 
commensurate solid and incommensurate solid, re- 
spectively. The vertical dash-dotted line marks the 
temperature above which the commensurate phase 
does not appear. 
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Figure 15. The argon-argon radial distribution functions for the system with = -18.6, recorded at 
r = 57 K and different values of /jj^ (shown in the figure). The vertical dotted lines mark the locations of 
subsequent maxima in a perfectly ordered commensurate phase. 



mensurate phase. Unfortunately, we cannot propose any reasonable explanation for the underestimation 
of the commensurate phase stabiUty by computer simulation. One can speculate that our model based 
on Lennard-Jones potential and standard mixing rules overestimates the tendency towards demtxing in 
submonolayer films. The commensurate phase is mixed, while the incommensurate xenon-like phase is 
demtxed. It is also possible, however, that x-ray diffraction data overestimate the range of xenon mole 
fractions corresponding to the commensurate phase. Note that within the region of coexisting commen- 
surate (C) and xenon-like incommensurate (IXxe) phases the paches of incommensurate phase may be 
quite small and hence escape detection. We recall that Villain and Moreira (s^ have also questioned the 
reliability of experimental results given in reference fisi] and suggested that the results were affected by 
metastability effects. 
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Figure 16. A comparison of the phase diagrams of Ar-Xe submonolayer films at two temperatures (given 
in the figure) resulting from the present Monte Carlo simulation and from the experiment 123]. 
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/\o (|)a30BoY noBeAiHKi/i SMimaHi/ix cy6MOHOLijapoBi/ix 
n;iiBOKAr-Xe Ha rpa(|)iTi 

A. naipi/iKeeB 

Biflfli/i M0fle;iK)BaHH5i (()i3WK0-xiMi4Hnx npoL^ecis, yHisepcuTeT Mapii Kropi-CK/iOflOBCbKoT, 
20031 J]Ki6n'\H, no;ibU(a 

BnKopMCTOByK)4n MeTOflw KOMn'iOTepHoro MOfle;iK)BaHHfl MeroflOM MoHie Kapyio y KaHOHi^HOMy i Be/iwKOMy 
KaHOHiHHOMy aHcaM6;iflx, mh o6roBopK)eMO n/iaB/ieHHfl i c|)opMyBaHHfl Bnop^iflKOBaHnx crpyKiyp BMiuiaHwx 
Ar-Xe cy6MOHOujapoBwx n/iiBOK Ha rpac|)iTi. 064wc;ieHH5i BUKOHyiOTbCfl 3 BWKOpwcTaHHflM flBO- i rpwBUMipHnx 
MOfleyibHux ci/iCTew. HoKasaHO, mo no3an;iou(MHHwR pyx He Bn/inBae na B/iacTHBOcri aflcop6oBaHOi n/iiBKn flo 
TUX nip, noKi/i 3ara;ibHa rycTuna ne crae 6;in3bK0K) flo MOHOuiapoBoro 3aBepujeHH5i. 3 iHuioro 6oKy, 6fli/i3bK0 
flo MOHOLuapoBoro 3aBepLijeHHfl, npocyBaHHfl nacTHHOK flo flpyroro mapy 3Ha4H0K) MipoK) Bn/iwBae na B;ia- 
CTMBOCTi 3MiLijaHnx H/iiBOK. noKa3aHO, iu,o cyMiui HOBHiCTK) 3MiujyeTbra B piflKiR c|)a3i i 3aMopo>KyeTbra y TBepfli 
c()a3M 3i CTpyKTypoK), mo 3a;iexnTb Bifl CK/iafly n/iiBKH. fl,n!f cy5M0H0Ujap0Bnx rycTHH, leMneparypa n/iaB/ieHHfl 
3MiHK)eTbra HeMOHOTOHHO 3i 3MiH0K) CK/iafly n/iiBKM. SoKpeMa, Tewneparypa n/iaB/ieHHn cnonaTKy spocrae 3 
pocTOM K0Hi4eHTpaL^iT KceHOHy 6;in3bK0 20%, noriM 3MeHUjyeTbCfl i flocarae MiniMyMy p^!f K0HL4eHTpaL4M kcbho- 
Hy 6;in3bK0 40%. fl/ia Bumwx K0HL4eHTpaL4iM KcenoHy T04Ka n/iaB/ieHHfl nocTynoBO 3pocTae flO leMneparyp, mo 
BiflnoBiflaiOTb n/iiBL^i 4HCTOro KceHOHy. TaKOx noKa3aHO, mo Tonoyiori^i i^asosm fliarpaM 3MiujaHwx n/iiBOK e 
HyT/inBOKD flO CK/iafly aflCop6oBaHwx LuapiB. 

K/iiOMOBi c;iOBa: afl,cop6u,m cyMimed, <pa30Bi nepexofli^, KOMn'tOTepne MOfle/iiOBaHHfi, nnaBJienHfi 
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